Quantum phases of interacting phonons in ion traps 
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The vibrations of a chain of trapped ions can be considered, under suitable experimental con- 
ditions, as an ensemble of interacting phonons, whose quantum dynamics is governed by a Bose- 
Hubbard Hamiltonian. In this work we study the quantum phases which appear in this system, and 
show that thermodynamical properties, such as critical parameters and critical exponents, can be 
measured in experiments with a limited number of ions. Besides that, interacting phonons in trapped 
ions offer us the possibility to access regimes which are difficult to study with ultracold bosons in 
optical lattices, like models with attractive or site-dependent phonon-phonon interactions. 



I. INTRODUCTION 



The interplay between atomic and many-body physics 
has proved to be an exciting research field in the last 
years. Experiments in atomic physics offer us the pos- 
sibility to find experimental realizations of theoretical 
models which were first proposed in the context of solid 
state physics. Many of these models are the key to under- 
standing a variety of phenomena in real materials such as 
quantum magnetism or superconductivity. For example 
cold bosonic atoms in optical lattices are a realization of 
the Bose-Hubbard Model (BHM) in a clean experimen- 
tal setup, where one can tune the value of interactions 
and observe quantum phase transitions in a controlled 
way [H . The main handicap of experiments with optical 
lattices is the fact that atoms are separated by optical 
wavelengths, and thus single particle addressing with op- 
tical means is severely limited by diffraction effects. 

Trapped ions are also an experimental system with po- 
tential applications to the quantum simulation of many- 
body problems @, 0, i, S, [| . It has the advantage that 
internal electronic or vibrational quantum states can be 
measured at the single particle level 0, 0] , since the dis- 
tance between ions is large enough to address them indi- 
vidually by optical means. In particular, we have recently 
shown that the vibrational modes of a chain of trapped 
ions under suitable experimental conditions follow the 
quantum dynamics of a BHM [3(. The interaction be- 
tween phonons is induced by the anharmonicities of an 
optical potential, which can be created by an off-resonant 
standing-wave. 

In this work we present a theoretical study of interact- 
ing phonons in trapped ions, and show the following re- 
sults: (i) The quantum phase transition between a super- 
fluid and a Mott insulator phase can be induced and ob- 
served in this system, (ii) Even though finite size-effects 
are important due to the finite length of the ion chain, 
properties corresponding to the thermodynamical limit 



can be accessed in experiments with a limited number 
of ions. These include critical exponents of correlation 
functions, and critical values of parameters in the Hamil- 
tonian. (iii) The ability to control phonon-phonon inter- 
actions allows us to realize models of interacting bosons 
which are difficult to reproduce in other experimental se- 
tups, like for example, BHM's with negative interactions, 
as well as models with site dependent interactions. 

The structure of the paper is the following. In sec- 
tion [TT] we derive the Bose-Hubbard model for phonons 
in a chain of trapped ions, in the presence of the anhar- 
monicities induced by an optical dipole potential. The 
DMRG algorithm that we have used to study numeri- 
cally this problem is summarized in section IIII1 In sec- 
tions IIVI and [Vl we study the quantum phases which 
correspond to repulsive and attractive phonon-phonon 
interactions, respectively. Section IVII is devoted to the 
case of a Bose-Hubbard model with site-dependent in- 
teractions. Finally in the last section we summarize our 
results and conclusions. 
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II. BOSE-HUBBARD MODEL OF PHONONS IN 
ION TRAPS 



In this section we show that under certain experimen- 
tal conditions, the dynamics of the vibrational modes of 
a chain of ions satisfies the Bose-Hubbard model of inter- 
acting phonons in a lattice Phonon number conser- 
vation is ensured whenever vibrational energies are much 
higher than other energy scales in the system. In this 
limit, physical processes which involve creation or de- 
struction of a phonon do not conserve particle number 
and are suppressed in much the same way as processes 
which do not conserve the number of electrons or atoms 
in low-energy physics. 
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A. Hamiltonian in the harmonic and phonon 
conserving approximation 

Let us start by writing the Hamiltonian that describes 
a chain of ions in a linear trap: 
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(1) 



N is the number of ions, and m is their mass. Pj and 
Rj are the momenta and the absolute positions of the 
ions, respectively. Vt is the trapping potential, which 
determines the ions' equilibrium positions. In this work 
we deal with two different situations. On one hand we 
consider the case of ions in a linear Paul trap, where they 
are confined by an overall trapping potential: 
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where u> a are the trapping frequencies in each spatial 
direction. On the other hand, we consider ions confined 
by an array of separate microtraps: 
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where Rj, a are the centers of each microtrap. Note that 
in ([3]) we assume that the confinement is strong enough, 
such that each ion feels only a single microtrap. 

The equilibrium positions of the ions are given by the 
minima of the trapping potential plus the Coulomb repul- 
sion. From now on, we choose the condition lu z 3> w x , u> y , 
such that the ion chain is along the z axis, with equi- 
librium positions given by z®. In the case of a linear 
Paul trap, described by Eq. ([2|), the equilibrium positions 
of the ions have to be calculated numerically. The dis- 
tance between ions is smaller at the center of the chain. 
In the case of independent microtraps (Eq. (J3|)), one 
can approximate the equilibrium positions by assuming 
that they correspond to the center of each microtrap: 
zf = Ri. z . This fact has strong implications for the 
phonon quantum dynamics, as we will see below. 

In the harmonic approximation, Hq is expanded up to 
second order in the displacements of the ions around the 
equilibrium positions, and we get a set of independent vi- 
brational modes corresponding to each spatial direction. 
The phonon number is a conserved quantity if the vibra- 
tional energies are the largest energy scale in the system. 
This condition can be met cither in the case of ions in 
individual microtraps, or in the case of the radial vibra- 
tions of ions in a linear trap, because the corresponding 
trapping frequencies can be increased without destroying 
the stability of the ion chain. For concreteness we restrict 
from now on to the case of vibrations in one of the radial 
directions, say x, but keep in mind that our results can 



be applied also to the axial vibrational modes if ions are 
in individual microtraps. 

The Hamiltonian that governs the dynamics of the ra- 
dial coordinates in the harmonic approximation reads: 
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where the displacements of the ions around the 

equilibrium positions, that is, simply Xi — Ri. x , and Pi tX 
the corresponding momenta. The second quantized form 
of this Hamiltonian is (we consider units such that h = 1): 
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(•5) 



a\ (ai) are creation (anihilation) operators for phonons 
in the radial direction. Harmonic corrections induced by 
the Coulomb interaction determine the effective trapping 
frequency, u> x ^ which depends on the ions' positions, as 
well as the tunneling amplitudes ti^\ 



UJ, 
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Eq. ([6]) yields an important result on the properties of 
phonons in trapped ions: the corrections to the local 
trapping energy, cu x i may depend on the position of the 
ions, in case the distance between ions changes along the 
chain. In [3| we have shown that in a linear ion trap, 
ions arrange themselves in a Coulomb chain, such that 
uj x _i is an effective harmonic confining potential for the 
phonons. On the contrary, in the case of an array of ion 
microtraps, the distances between ions can be considered 
to be approximately constant, and thus this confining ef- 
fect does not take place. 

Before going any further, let us study under which con- 
ditions phonon nonconserving terms of the form (a^a,,-, 
ajat), can be neglected in Eq. ([5]). We define the param- 
eter: 



P x = e7(m^^), 



(8) 



where do is the distance between ions. Since we will be 
interested in the limit f3 x -C 1, we choose do to be the 
minimum distance between ions in the case of a linear 
Paul trap. Phonon tunneling terms Uj are of the order 
of t, defined by: 



t = f3 x uj x /2. 



(9) 
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Since phonon nonconserving terms rotate fast in ([5]), we 
can neglect them in a rotating wave approximation if: 

t/u x = x /2 « 1, (10) 

such that Eq. ([5]) takes the form of tight binding Hamil- 
tonian with hopping terms tij(a\aj + h.c). 

In our numerical calculations, we will parametrize the 
tunneling of phonons between sites by the parameter t, 
which corresponds, due to the definition of /3 X , to the 
highest value of the tunneling along a chain in the case of 
a linear Paul trap, and to the tunneling between nearest- 
neighbors in the case of an array of microtraps. 

B. Phonon— phonon interactions 

Anharmonic terms in the vibrational Hamiltonian are 
interpreted as phonon-phonon interactions, and can be 
induced by placing the ions at the minimum or maximum 
of the optical dipole potential created by an off-resonant 
standing wave along x: 

N 

H sw = Fj2cos 2 (kx i + ^6). (11) 
i=l 

F is the amplitude of the dipole potential, and 5 deter- 
mines the position of the ions relative to the standing 
wave. We define the Lamb-Dicke parameter r\ — kxo, 
where k is the wave-vector of the standing wave lasers, 
and xq is the ground-state size of the radial trapping 
potential. The only relevant cases for us are 5 = (max- 
imum of the optical potential), and 5=1 (minimum). 
Under the condition r\ -C 1, we can write H sw as a se- 
ries around x% = 0. The term that is quadratic in xi 
in Eq. (jlip can be included in the harmonic vibrational 
Hamiltonian just by redefining the global radial trapping 
frequency: 

ul - lu x (lj x - {-fAifF) . (12) 

In the case 5 = 0, condition rj 2 F <C lo x has to be fulfilled, 
such that the radial trapping frequency is not strongly 
suppressed by the standing wave, and the system remains 
in the phonon number conserving regime. Under this 
condition the only relevant term is thus the quartic one: 

We can neglect nonconserving phonon terms again under 
the condition _Fr/ 4 CWj. In this way we get, finally, the 
promised BHM for phonons: 

JV 

# X BHM = X ■'>■<■■> 

i>j 

N N 

+ "^2(LO x +LO Xtl )alai + Uj2al 2 a 2 . (14) 

i=l i=l 



The on-site interaction is given by: 

U = 2(~1) S F V 4 . (15) 

Thus, it is repulsive or attractive depending on whether 
the ions are placed at a minimum or maximum of the 
standing-wave, respectively. 

III. NUMERICAL METHOD 

The Bose-Hubbard Model with tunneling between 
nearest-neighbors has been thoroughly studied in the 
past [1, Hof . It has recently received considerable at- 
tention because it describes experiments with ultracold 
atoms in optical lattices. In general, we expect the same 
phenomenology to appear in our problem, such as, for ex- 
ample, a superfluid-Mott insulator quantum phase tran- 
sition. However, the situation of phonons in ion traps 
presents a few peculiarities that deserve a careful anal- 
ysis: the effects of long-range tunneling in (jT4j) , finite 
size effects, as well as the possibility of having attractive 
interactions. 

To handle this many-body problem numerically we 
use the Density Matrix Renormalization Group method 
(DMRG) [ll|, which has proved to be a quasi-exact 
method in quantum chains. In particular we use the 
finite-size algorithm for open boundary conditions. Our 
problem is defined in the microcanonical ensemble, that 
is, we find the minimum energy state within the Hilbcrt 
subspace with a given number of total phonons N p ^. For 
this reason, we have implemented a DMRG code which 
uses the total phonon number as a good quantum num- 
ber and projects the problem into the corresponding sub- 
space at each step in the algorithm [l2T |. 

To keep a finite dimensional Hibert space, we truncate 
the number of phonons in each ion, and define a max- 
imum value n max , which is usually taken to be of the 
order of 6(n), with (n) the mean phonon number. The 
number of eigenstates of the reduced density matrix that 
we keep at each step is always in the range 80 — 100. 

We have checked the accuracy of our method by com- 
paring our numerical calculations with the exact solu- 
tion in the case of a system of non-interacting (U = 0) 
phonons, where the ground state at zero temperature is 
a condensate of phonons in the lowest energy vibrational 
state. We have also compared our numerical results with 
exact diagonalizations of Eq. (| 14[) with up to N = 5 ions. 
In both cases we found agreement between DMRG and 
the exact results up to machine accuracy 5E ~ 10~ 14 . 

The relevant experimental parameters of our phonon- 
Hubbard model are discussed in the Ref. Q. Typically 
we could choose the minimum distance between ions do = 
5 fim and f3 x = 2 x 10~ 2 <C 1. Then we have uj z « 
177 kHz and oj x » 12.5 MHz for a string of ions with 
N = 50. The number of phonons is N p h = N for the 
superfluid and Mott-insulator phases and N p h = N/2 for 
the Tonks-gas phase. In the end we discuss N p h = 2N 
for a special case with site-dependent interactions. All 
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our calculations are for the ground state, i.e. at zero 
temperature. 

Following the discussion below Eqs. ([6j [7]), one ex- 
pects to find significant differences between the cases of 
phonons in ions trapped in a linear trap (Coulomb chain), 
and phonons in an array of ion microtraps. Finite size 
and inhomogeneity effects are indeed much more impor- 
tant in the linear trap case, since harmonic Coulomb cor- 
rections induce an effective harmonic trapping for the 
phonon field. For this reason, we always study these two 
cases separately, in the different quantum phases that we 
will explore in what follows. 




FIG. 1: (a) Mean phonon number (rtj), and (b) fluctuations 
Srij , at each ion in the phonon ground state in an array of mi- 
crotraps. Number of ions N = 50, and total phonon number 
iVph = 50. 



IV. BOSE-HUBBARD MODEL WITH 
REPULSIVE INTERACTIONS: U > 

We study first the quantum phases of phonons with 
U > 0, and both commensurate and incommensurate 
total phonon number. In this section we present results 
for a chain with N = 50 ions, and total phonon number 
N p h — N in the commensurate case, or N p ^ = N/2 in 
the incommensurate case. 

The local observables that we consider are the num- 
ber of phonons at each site, rij 

fluctuations 



[a]aj) 



Srij = 



as well as its 
2 . Two- point correla- 



V V3 

tion functions have to be defined carefully, to take into 
account finite size effects, for example, the variations of 
the density of phonons along the chain. Correlations in 
the number of phonons are given by: 



(16) 



A suitable definition of correlations that are non- 
diagonal in the phonon number basis is the following one 



13]: 



( a i a j) 

\/{ni){nj) 



(17) 



such that correlations are rescaled by local values of the 
phonon density. The rescaling is inspired by the decom- 
position of the phonon field in density and phase opera- 
tors, a,j = ^Jrije'^i , which is the starting point for the 
Luttinger th eory of the weakly interacting bosonic super- 
fluid 0, GJ, G]f . 

In Figs. [T] and [2l we plot the local density and its 
fluctuations in the cases of an array of microtraps, and a 
linear Paul trap, respectively. These figures show a sig- 
nature of the different phases which can be observed in 
the model defined by the Hamiltonian (fl4)) . Figs. []Ja) 
and[^a), in particular, show the variation of the den- 
sity of phonons along the chain. The evolution of the 
density profile shows the transition from the phonon su- 
perfluid to the Mott-insulating phase. When t ^S> U, the 
ground state of the system is a condensate such that all 
the phonons occupy the lowest energy vibrational mode. 
In the case of a linear Paul trap, phonons are confined 
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FIG. 2: (a) Mean phonon number and (b) fluctuations, at 
each ion in the phonon ground state in a linear Paul trap. 
N = 50, Ap h = 50. 



in the center of the chain, due to the effective trapping 
potential induced by the nonconstant ion-ion distance. 
At U ^> t, the ground state is a phonon Mott insulator 
with approximately one phonon per site, and no phonon 
number fluctuations. Note that due to the effective har- 
monic trapping potential, the Mott phase in the whole 
chain is reached for lower values of U in the case of the 
array of microtraps (Fig. [1} than in the linear Paul trap 
case (Fig. [2j. 

In the following subsections, we study these two quan- 
tum phases separately, paying particular attention to 
their correlation functions. 



A. Superfluid phase 

When the tunneling dominates the on-site interac- 
tions, the system is in the superfluid phase [13] ■ The 
non-interacting ground state is given by a condensate 
solution in which the A p h phonons are in the lowest vi- 
brational mode: 



phi 



5> 



|0), (18) 



where is the wave-function of the lowest energy vi- 
brational mode. Interactions suppress long range order 
in ID, even in the weak interacting limit, U <C t, in which 
Luttinger liquid theory allows us to make predictions on 
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FIG. 3: Correlation functions (a) Cf^ and (b) C"™ as a 
function of coordinate j, at the superfiuid phase (U/t = 0.1) 
in an array of ion microtraps. We choose io = 26 (center 
of the chain), N = N p h = 50. The dotted and solid lines are 
numerical results and fittings in the region where the functions 
show algebraic decay. 



the scaling of correlation functions: 

C%<x\i-j\- a , 

C^cx|<-j|- a . (19) 
where a depends on the parameters of the model: 



a oc 



' U/t 
n 



(20) 



In deriving (|19[) one has to neglect phonon tunneling be- 
yond nearest-neighbor ions, and assume an homogeneous 
system LJ]. In the following we will check if Luttinger 



theory describes also our numerical results in the case of 
phonons in a chain of trapped ions, by fitting our results 
to the form ifTT 



(1) Array of microtraps. We start with the case of 
the superfiuid phase in an array of ion microtraps, see 
Fig. [3l Correlation functions Cf°j and C™" decay alge- 
braically in an intermediate range of ion-ion separations, 
with exponents which satisfy the predictions of Luttinger 
theory. In particular, the evolution of a in Eq. (|19[) is 
well described by the Luttinger liquid scaling law ([20]) . 
as shown in Fig. Q] 

(2) Linear ion trap. In the case of ions in a linear 
Paul trap, finite size effects play a more important role, 
because of the inhomogeneities of the on-site phonon en- 
ergy. Correlation functions still decay algebraically for 
short distances in the superfiuid regime, but boundary 
effects spoil this behavior at large separations between 
ions. In the algebraic regime, exponents are close to those 
predicted by Luttinger theory in the homogeneous case, 
see Figs. [5]and|6l 

Due to the localization of phonons as we increase U/t, 
a Mott insulator phase appears first at the sides of the 
ions chain, which coexists with a superfiuid core at the 
center. This coexistence of the phases can be observed in 
the correlation functions, which show regions of algebraic 
or exponential decay, as shown in Fig. [7J 
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FIG. 4: The evolution of the exponent a of C™" • as a function 
of U /t in the regime of superfiuid phase, in an array of ion 
microtraps (io = 26, N = iV p h = 50). The dotted line is from 
numerical data, and the solid line the fitting result from the 

Luttinger liquid, that is, a m A^J^-, where the coefficient 
Aw 1.68. 





FIG. 5: 



Correlation functions (a) C^j 



and (b) at the 

superfiuid phase U/t — 0.2 in a linear Paul trap (N = N p h = 
50). The dotted lines are numerical data, and the solid lines 
are fittings in the region where the correlations decay alge- 
braically. 
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FIG. 6: The evolution of the parameters a in Ci//j with U/t. 
The dotted line is the values from the correlation, and the 



solid line is the fitting data from a 
ing to the Luttinger liquid theory. 
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FIG. 7: Correlation functions (a) Cf^j 

linear Paul trap when U/t = 2 (AT = JV ph = 50). Both of 
them show the coexistence of superfluid and Mott-insulator 
phases. In (a) and (b), the insets (i) show the exponential 
decay in the region of the Mott phase, and the insets (ii) 
show the occupation number and fluctuations at the same 
parameters. The exponents a of the algebraic decay are also 
given in the figures. The dotted and solid lines are numerical 
and fitting data, respectively. 



B. Mott-insulator phase 

In the commensurate case, by increasing the on-site 
interaction U, a quantum phase transition from a su- 
perfluid to a Mott-insulator state takes place at about 
U ~ 2t. In the limit in which interaction dominates over 
hopping, the ground state for a commensurate filling of n 
particles per site is simply a product state of local phonon 
Fock states, 



MI I 



N 1 

n4y(4) B io>. 



(21) 



In the Mott insulator phase correlations decay exponen- 
tially with distance, C, 



oc e"' 1 ^ 1 ''^, where £ is the 
correlation length. The correlation length diverges when 
approaching the quantum phase transition. 

(1) Array of ion microtraps. In Fig. [5] (a) we plot 
the correlation functions in the phonon Mott phase in 
the case of an array of ion microtraps. These curves can 



FIG. 8: (a) Correlation functions Cf °,- at U/t ■ 



■ 3 and CI 



at U/t = 2.4 decay exponentially with the coordinate j at the 
Mott-insulator phase for an array of microtraps (N = iV p h = 
50, io — 26). (b) The inverse of the correlation lengths in 
Cf^j and C™y The circle markers represent Cf^j and the 
square markers CfJ^j. The solid lines are the fitting data. 
The extrapolation in (b) show that the critical point in the 
microtraps is U c /t « 1.55. 



FIG. 9: The exponential decay correlations and C™j at 
the Mott phase U/t — 2.8 for a linear Paul trap (N = N p ^ = 
50, io = 26). 



be fitted to an exponential decay, and the correspond 
correlation lengths are plotted in Fig. [5] (b) as a func- 
tion of the interaction strength. Due to the finite size 
of the system, £ does not diverge at the critical value of 
U. However, the extrapolation of the curves in the linear 
regime allows us to estimate the critical point, which lies 
at U c /t w 1.55. This critical value is smaller than the one 
in the BHM with tunneling between nearest -neighbors 
only, Ul m /t « 2. The condition U c < U™ is due to the 
frustration induced by hopping between next-nearest- 
neighbors, which makes the superfluid phase more un- 
stable against the effect of on-site interactions. 

(2) Linear ion trap. In the case of a linear trap, the 
behavior of spatial correlations is similar, see Fig. [9] C"™ 
is difficult to fit due to the few points with exponential 
decay, therefore we only plot the correlation length cor- 
responding to Cfj . Due to the effective phonon trapping 
potential, the Mott insulator and superfluid phases co- 
exist in a range of values of U (Fig. [7|). For this reason, 
in the case of a linear Paul trap, one cannot follow the 
extrapolation procedure of Fig. [5] to find a critical value 
of the interaction. 

Finally, in the Mott-insulator phase the long-range 



hopping terms which decay like l/\z° 



.013 



play a major 



role, since they induce a peculiar long-range correlation in 
this phase. In Fig. [TU1 we show that indeed also be- 
haves like l/|Zj° — z*j\ 3 at long distances. The existence of 
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FIG. 10: The correlation C?^ at the Mott phase U/t = 6 in 
an array of microtraps (circle markers) and a linear Paul trap 
(square markers), respectively. N = iV p h = 50, io = 26. The 
inset shows the power-law decay with exponent a ~ 3 for the 
microtraps only. 

power-law decay in correlation functions of non-critical 
systems due to long-range interactions was also observed 
in the case of spin models in trapped ions, see the dis- 
cussion in Ref. Hi. 



C. Tonks-gas phase 

We turn now to the incommensurate filling case, where 
in the limit U ^> t the system forms a Tonks-Girardeau 
gas, which can be described in terms of effective free 
fcrmions. A Tonks-Girardeau gas has recently been real- 
ized in an experiment with ultracold bosons in an optical 
lattice QJ]. 

(1) Array of ion microtraps. We have studied numeri- 
cally the Tonks-Girardeau regime in the case of phonons 
in ion traps, starting with the case of an array of ion mi- 
crotraps with N = 50 sites and N p h = 25 phonons, that 
is, 1/2 filling. The density of phonons evolves from a 
superfluid to a Tonks-gas profile when increasing the in- 
teraction, and at the end it approaches a constant value 
of 1/2. Correlation functions decay algebraically, with 
an exponent that approaches a ~ 0.58 for large inter- 
actions (see Fig. HTj) . Note that a deviates from 1/2, 
which is the value that corresponds to a Tonks gas with 
nearest-neighbor tunneling only. The deviation can be 
explained by the mapping from the BHM model (|14p with 
U 3> t to an XY model with antiferromagnetic interac- 
tions of the form J^j = J/\i — j| 3 . The long-range terms 
in the antiferromagnetic interaction induces a change in 
the exponent of the correlation functions, as shown with 
the numerical calculations of our previous work on spin 
models in ion traps 

(2) Linear ion trap. We study now the case of phonons 
in a linear Paul trap under the same conditions, see Fig. 
H"2"l Correlation functions decay algebraically, with an 
exponent that is extrapolated to a = 0.53 in the limit of 




FIG. 11: (a) Densities of phonons in an array of microtraps 
with N = 50, N p h = 25. (b) Fluctuations at the same condi- 
tions, (c) The correlation function Cf^j (io = 26) at Tonks- 
gas phase U/t = 10 with exponent a « 0.54, where the dotted 
and solid lines are numerical and fitting data, respectively, 
(d) Evolution of exponent a of C^j with the interaction U/t, 
which would approach a ~ 0.58. 




FIG. 12: (a) and (b) show densities and fluctuations of 
phonons, respectively, in a linear Paul trap with N = 50 
and Nph = 50; (c) The correlation Cf^j at the Tonks-gas 
phase U/t = 6, whose exponent is a w 0.48; (d) Evolution of 
the exponents a of C^j, approaching a w 0.53. In (c) the 
dotted and solid lines represent numerical and fitting data, 
respectively. 



strong interactions. This result can also be explained by 
the mapping to the XY model, and coincides with the 
result found in 

In order to test if the system is really in the Tonks- 
gas phase, we introduce the quantity (O) — (J2i n i( n i ~ 
1))/N, which measures the probability of phonon oc- 
cupancies larger than one. In the Tonks-gas regime 
(O) ~ 0. The parameter (O) as a function of the in- 
teraction U is plotted in Fig. [T31 showing the continuous 
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FIG. 13: The evolution ol the parameter (O) with the inter- 
action U /t in an array of microtraps (solid line) and a linear 
Paul trap (dashed line). N = 50, N ph = 25. 



evolution into the Tonks-gas regime. 

Our results are consistent with the behaviors observed 
in optical lattices . The numerical analysis shows that 
phonons in ion traps are also a good candidate like atoms 
in optical lattices for studying Tonks gases. 



V. BOSE-HUBBARD MODEL WITH 
ATTRACTIVE INTERACTIONS: U < 

The BHM with attractive interactions in optical lat- 
tices has been the focus of recent theoretical studies . 
The sign of phonon-phonon interactions in trapped ions 
can be made negative simply by changing the relative po- 
sition of the standing-wave relative to the ion chain. For 
a qualitative understanding of this model, it is useful to 
consider a Bose-Einstein condensate in a double well po- 
tential H3| • In a symmetric potential in the absence 
of tunneling, energy is decreased when bosons accumu- 
late in one of the wells. When tunneling is switched on, 
the ground state of the system is a linear superposition 
of states with all the bosons placed in one of the wells, 
showing large phonon number fluctuations. 

The increase of phonon number fluctuations in our 
model when we switch on a negative interaction is also 
shown in our numerical calculations. In Fig. 1141 the 
density at the center of the ion chain and its fluctua- 
tions increase with the magnitude of the interaction for 
N = 10 ions and N p h — 10 phonons with open boundary 
conditions. Due to the open boundary condition and the 
symmetry of the potential, the phonons tend to collect 
themselves on one of the two sites at the center when 
increasing \U\ with an even number of sites. The ground 
state is then a superposition of N p h phonons on site N/2 
and Nph phonons on site N/ 2 + 1 . 

When |Z7| is large enough, our numerical calculations 
yield a ground state with all the phonons in a single ion, 
such that the spatial symmetry of the problem is broken. 
This effect is an artifact of the DMRG calculation, due 
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FIG. 14: (a) Density of phonons and (b) phonon number 
fluctuations in an array of microtraps with N = 10, N p h = 10, 
and negative on-site interactions . (c) Density of phonons and 
(d) phonon number fluctuations in a linear Paul trap under 
the same conditions. 




FIG. 15: Evolution of the order parameter (O) with the ratio 
U/t for an array of ion microtraps (solid line) and a linear 
Paul trap (dashed line). N = N ph = 10. 



to the small energy difference between the exact ground 
state of the system and the one which breaks the spatial 
symmetry. Thus, in order to study properly the phonon 
phases with negative interaction, it is convenient to de- 
fine the following order parameter, whose value is inde- 
pendent on the breaking of the spatial symmetry in the 
problem: 



(22) 



In Fig. Qj)] we show the evolution of this quantity, which 
shows a sudden increase for negative interactions. 
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FIG. 16: The phonon density (rij) and phonon number fluc- 
tuations, Srij, in an array of ion microtraps, with the on-site 
interactions defined by Eq. [f23]l. U/t = 40, TV = 50, and 
N ph =2N= 100. 



VI. BOSE-HUBBARD MODEL WITH 
SITE-DEPENDENT INTERACTIONS 

Phonons in trapped ions have a higher controllabil- 
ity than ultracold neutral atoms in optical lattices, due 
to the possibility of individual addressing. In particu- 
lar, on-site interactions can be induced in such a way 
that they depend on the ion position. In this section we 
present a model which shows how this possibility can be 
exploited for the engineering of quantum phases. 

Let us consider repulsive on-site interactions which 
vary over the ion chain in the following way: 

Ui = U dd = U, i odd, 

Ui = [/even = 2C7, i even, (23) 

We focus on the case with filling factor 2, A p h = 2N, 
in the regime where interactions dominate over tunnel- 
ing, U/t ^ 1. In the limit t — 0, the ground state of 
this model is highly degenerate. For instance, in a chain 
with two sites the ground state manifold in the Fock ba- 
sis spans the states |2, 2) and |3, 1). In a chain with even 
number N of sites and 7V p h = 2N, the ground state de- 

generacy is I ^ I . 

Our DMRG algorithm allows us to calculate the den- 
sity and fluctuations in the phonon number, which are 
shown in Fig. [TH for the case where the interactions 
defined by Eq. (|23[) are induced on an array of ion mi- 
crotraps. In the ground state of the chain, the number of 
phonons fluctuates between |2) and |3), and |2) and |1), 
in odd and even sites, respectively. 

This model can be understood in the hard-core boson 
limit by introducing a spin representation, which is valid 
near filling factor 2. At each site in the ion chain, we 
define a two level system by means of the following rule: 

|Q)i = |2) i) \l)i = \3)i, iodd, 

16)4 = 11),, \l)i = \2) h ieven. (24) 

Where |0) and |1) are the two levels which define the spin 
states of the spin representation of the hard-core bosons. 



FIG. 17: The correlation function (cr+Oj) defined by p4)l . 
io = 26, N = 50, A ph = 100, U/t = 40. ° 



Spin and phonon annihilation operators satisfy: 

<rf = iaj, i odd, 

o-+ = -^=aj, ieven. (25) 

where the equality is understood to hold within the 
ground state manifold. In terms of this operators, the 
Hamiltonian of the system is described by an XY model 
(for simplicity we consider here the nearest-neighbor 
case): 

H = jJ2( a t^ + i+h.c), (26) 

i 

with J — v6t- Under the condition N p h/N = 2, the 
ground state of our hard-core boson Hamiltonian corre- 
sponds to the solution of the XY model (|26[) with the 
constraint V\ a t = 0- 

Spin-spin correlation functions are related to the 
bosonic correlations functions of the BHM by the rela- 
tion Eq. (|2"B"|) . In Fig. [T7] we plot (crfcTj) calculated 
by means of correlation functions of hard-bosons. This 
correlation function shows an algebraic decay for short 
distances, which is spoiled for long separations between 
ions due to boundary effects. The exponent a ~ 0.56, 
differs from the one that we expect from the mapping 
to the XY model, that is, a = 0.5, due to the effect of 
further than nearest-neighbor interactions terms, which 
we have neglected. 

Another interesting possibility is to study the phase di- 
agram of the BHM with filling factor 2, and alternating 
interactions, beyond the XY point U even = 2U dd- If we 
change the ratio U even /U dd in the vicinity of this point, 
the ground state degeneracy with zero tunneling is lifted. 
The system is in a Mott insulator phase, with constant 
phonon density if U cvcn > 2U dd, |2, 2, 2, . . . ), or alternat- 
ing occupation numbers if U cvcn < 2U dd , |3, 1,3,1,...). 
The phase diagram as a function of the ratio U even /U dd, 
shows two gaped regions, separated by a single critical 
point, which corresponds to the XY limit studied above. 
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FIG. 18: The energy gap AE between the ground state and 
the first excited state as a function of U cvcn /U dd- U dd/t = 
40. For simplicity, we consider here a small system of ions 
TV = 6. 

This is shown in Fig. [IH where we calculate the energy 
gap AE for small chains (N = 6). 

To generalize, the same also takes place with other 
distributions of the on-site interactions, whenever the 
number of sites where phonons interact with U is the 
same as the number of sites with 2U interaction. For 
example, the chain can be divided in two regions, left 
and right, such that the interactions depend on the site 
in the way U, U, 2U, 2U, that is, Ui — f7i e ft = U and 
Ui = f/ r i gh t - 2U, 



VII. CONCLUSIONS 

In conclusion, we have studied the quantum phases of 
interacting phonons in ion traps. The superfluid-Mott 
insulator quantum phase transition can be detected by 
the evolution of the phonon density profile, as well as by 
the divergence of the correlation length near the quan- 
tum critical point. Although boundary effects are im- 
portant, specially in the case of a Coulomb chain of ions, 
correlation functions show a similar behavior as those 
of systems in the thermodynamical limit. For example, 
Luttingcr liquid theory gives an approximate description 
of the algebraic decay of correlations in the superfluid 
regime. 

We have also shown that the ability to control phonon- 
phonon interactions allows us to study a variety of sit- 
uations like attractive interactions, where a phase with 
large phonon number fluctuations takes place. The abil- 
ity to tune locally the value of the on-site interactions 
also leads to the realization of new exciting models, where 
the degeneracy of the classical ground state can be tuned 
by choosing properly the value of the phonon-phonon in- 
teractions. 

Work supported by the Deutscher Akademischer Aus- 
tausch Dienst, SCALA and CONQUEST projects and 
the DFG. 
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